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ANALYSIS  OF  THE  IIP  ICEBERG  DETERIORATION  MODEL 


ABSTRACT 

The  International  Ice  Patrol  has  developed  a  model  of  iceberg  deterioration 
that  plays  a  critical  role  in  the  overall  modeling  of  iceberg  locations  and 
determining  the  Limits  of  All  Known  Ice.  The  deterioration  model,  along 
with  an  iceberg  drift  model,  is  used  to  determine  the  location  and  state  of 
any  reported  iceberg  until  it  is  resighted.  The  deterioration  model  depends 
on  the  reported  size  of  the  iceberg  and  a  number  of  environmental 
parameters.  This  report  examines  the  structure  of  the  deterioration  model, 
reviews  empirical  evaluations  of  the  model,  and  conducts  an  analytical  and 
experimental  sensitivity  analysis  of  the  model  outputs  with  respect  to  the 
model  parameters.  The  result  of  this  analysis  suggests  that  the 
deterioration  model  is  a  reasonable  representation  of  the  actual 
deterioration  process,  lacking  only  a  good  mechanism  for  modeling  the 
calving  phenomenon.  The  sensitivity  analysis  suggests  tlmt  there  is  little 
need  for  further  accuracy  in  obtaining  environmental  parameters.  The 
single  most  important  factor  in  the  application  of  the  model  is  the  initial 
estimate  of  iceberg  size  and  the  corresponding  waterline  length. 


INTRODUCTIpN 

Objective. 

The  IIP  uses  an  Iceberg  Deterioration  Model  to  emulate  the  approximate 
deterioration  behavior  of  icebergs  in  the  presence  of  known  environmental  conditions. 
The  purpose  of  this  report  is  to  review  the  structure  of  the  existing  model  and  to  examine 
the  sensitivity  of  its  output  with  respect  to  changes  in  input  parameters.  Tliis  analysis 
identifies  which  input  parameters  require  the  most  attention  with  respect  to  accuracy  of 
their  estimates  and  identifies  areas  where  potential  model  enliancements  are  appropriate. 


Background. 


The  iceberg  deterioration  model  used  by  IIP,  based  on  the  work  of  Wliite, 
Spaulding,  and  Gominlio  (1980),  was  completed  and  initially  tested  in  1983  (Anderson, 
1983).  The  model  considers  four  forms  of  deterioration:  insolation  (sun  heating),  buoyant 
convection  (vertical  circulation  of  the  water),  wind  forced  convection  (drift  movement 
through  the  water),  and  wave  induced  deterioration  (wave  wasliing  of  the  subaerial 
surface).  Four  equations  determine  the  melt  due  to  each  of  these  processes  wliich  are 


additive.  Input  data  include:  iceberg  position;  iceberg  size;  sea  surface  temperature;  wave 
height;  and  wave  period.  Progressive  deterioration  of  the  iceberg  is  quantified  by  its 
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waterline  “length”  by  definition.  There  are  four  size  categories:  growler,  small,  medium, 
and  large  with  no  upper  limit.  The  model  used  by  IIP  does  not  include  deterioration  due 
to  calving  of  overhanging  ice  slabs  (Anderson,  1983).  Deterioration  due  to  calving 
depends  on  the  thickness  of  the  overhanging  slab  which  is  not  generally  available.  Not 
including  the  deterioration  due  to  calving  underestimates  the  degree  of  deterioration  in  a 
given  time  period.  The  effect  of  the  four  forms  of  deterioration  are  illustrated  in  Figure  1 


Modelled  Iceberg  Deterioration  Processes.  This  figure  depicts  four  processes  used  by 
the  IIP  deterioration  model  to  "melt''  icebergs.  The  four  processes,  which  are  labelled  in  order  of 
importance  are:  (1)  wave  erosion;  (2)  forced  heat  convection;  (3)  buoyant  heat  convection;  and  (4) 
solar  radiation. 

Figure  1.  Modeled  Iceberg  Deterioration  Process  (Hanson,  1987). 
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Iceberg  Deterioration  Model. 

The  Iceberg  Deterioration  Model  (Anderson,  1983)  consists  of  four  additive 
components  that  contribute  to  deterioration  in  waterline  length  (cm/day)  by  way  of  heat 
transfer:  j’j  represents  deterioration  due  to  insolation;  >^2  represents  deterioration  due  to 
buoyant  convection;  represents  deterioration  due  to  wind  forced  convection;  and  y4 
represents  wave  induced  deterioration.  Seven  parameters  (related  to  the  notation  of 
Anderson  (1983))  are  used  to  characterize  these  four  components: 

x^  =  SST  (Sea  Surface  Temperature  in  "C) 

Xj  =  XAMP  (Wave  Height  in  cm) 

Xj  =  RELSPD  (Relative  Speed  in  cm/s) 

X4  =  RLEN  (Waterline  length  in  m) 

Xj  =  IPER  (Wave  Period  in  s). 

Xg  =  ZTIME  (in  halfdays) 

X7  =  WEATHER  (1  =  CLOUDY,  2  =  CLEAR). 

The  weather  in  the  IIP  operating  area  is  seldom  clear,  therefore  x^  =  1.  To  obtain 
melt  rate  per  day,  Xg  =  2.  The  resulting  model  for  the  total  melt  rate  Y  in  cm/day  is  given 
by  equation  (1) 

Y=yi+y2+y3+y4 


where 


Insolation: 

jj  =  X6X7/100  (2) 

=  .02 

Buoyant  convection 

J2  =  0.274*(2.78T  +  0.47TV0.5*Xj/100 

(T  =  X,  -  Iceberg  Surface  Temperature) 

(Iceberg  Surface  Temperature  was  held  constant  at  -\°C) 

=  0.0012878x,^  +0.0101928X,  +  0.008905  (3) 

Wind forced  convection 

J3=FC  (x,  +  1)*0.5*x,/100 

=  FC(x,  +  1)/100  (4) 

where 

FC  =  (0.934-0.202log,oX4  )x3  (x,  <  25) 

FC  =  (0.66-0. 15 Hog, 0X4 )(x3-25)+25(0.934-0.2021og, 0X4)  (X3  >  25) 
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and 


Wave  induced  deterioration 

*  =  0.000146x2  (2/x2)02*24*3600*[(x,  +  1)/X5]*0.5*X6/100 

=  0.144901x/(Xi  +  1)/x5  (5) 

The  model  equations  applied  with  a  6'  wave  height,  10  sec  wave  period,  and  25 
cm/sec  relative  velocity  (SST  unknown)  result  in  84  percent  of  the  iceberg’s  deterioration 
being  attributable  to  wave  induced  melting;  approximately  14  percent  due  to  wind  drift  of 
the  iceberg  relative  to  the  water  bathing  it;  and  less  than  2  percent  related  to  sun  heating 
and  vertical  circulation  along  the  iceberg's  submerged  surface. 

In  the  application  of  the  model,  the  maximum  waterline  length  (x^)  for  the  size 
category  reported  is  used  initially.  If  an  iceberg  is  resighted,  the  waterline  length  is  set  to 
the  maximum  for  the  size  reported  in  the  resighting  eliminating  any  deterioration  which 
may  have  occurred  if  the  resighting  indicates  the  same  or  larger  size.  Icebergs  for  wliich 
no  size  category  is  reported  are  assumed  to  be  medium  icebergs.  Note  that  the  only 
melting  component  that  depends  on  the  size  of  the  iceberg  (waterline  length)  is  wind 
forced  convection  melting.  Thus,  deterioration  can  be  computed  after  the  waterline  length 
has  decreased  to  zero.  Icebergs  are  removed  fi'om  the  model  when  125  percent  of  their 
original  waterline  length  has  been  melted  if  they  remain  witliin  the  bounds  of  all  known 
ice,  unless  they  are  limit  setting  icebergs  for  the  region  of  all  known  ice  in  which  case  they 
are  retained  by  the  deterioration  model  until  150  percent  of  their  waterline  length  has  been 
melted. 


QUALITATIVE  AND  EMPIRICAL  EVALUATION  OF  THE  MODEL 

Anderson  (1983)  conducted  an  initial  mathematical  analysis  of  the  model:  with  all 
other  parameters  held  constant,  a  100  meter  length  iceberg  took  179  days  to  melt  in  -1°C 
water  compared  with  20.5  days  in  3°C  water.  Anderson  concluded  that  input  data  errors 
of  1°C  variance  firom  actual  sea  surface  temperatures  in  tliis  temperature  range  can 
produce  melt  errors  on  the  order  of  40  days  for  this  100  meter  berg  length.  These  results 
suggest  that  the  sea  surface  temperature  is  the  most  critical  parameter.  However,  other 
parameters  were  not  examined  and  the  results  only  hold  in  the  range  of  SST  considered  in 
the  analysis. 

Several  studies  by  Venkatesh,  El-Tahan,  and  Mitten  (1985),  Venkatesh  (1986), 
and  El-Tahan,  Venkatesh,  and  El-Tahan  (1987)  compared  model  performance  with 
observed  deterioration  of  several  icebergs  using  observed  oceanograpliic  and 
meteorological  data.  They  also  developed  refined  size  estimates.  Using  these  data,  they 
found  good  agreement  between  the  model  results  and  actual  deterioration. 

In  1987  the  IIP  conducted  a  deterioration  study  using  6  icebergs,  observed  and 
tracked  by  a  surface  vessel  (Hanson,  1987).  The  time  of  observation  on  each  iceberg 
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ranged  from  2.1  days  to  6.3  days.  The  objectives  of  this  study  were  to  compare  iceberg 
deterioration  predictions  derived  from  observed  environmental  to  predictions  using  system 
(FNMOC)  data.  FNMOC  SST  data  were  an  average  of  1.3°C  colder  than  that  actually 
observed,  the  wave  heights  averaged  0.9  meters  higher  than  observed,  and  the  periods 
were  on  average  4.6  seconds  greater  than  observed.  The  6  iceberg  cluster  averaged  379 
cm/day  melt  rate  of  their  waterline  length  using  observed  wave  erosion  values,  while  using 
the  system  operational  data  provided  by  FNMOC  produced  a  melt  rate  on  average  of  531 
cm/day.  The  overestimation  of  the  predicted  wave  height  was  identified  as  the  primary 
cause  of  the  significant  overestiniation  of  the  melt  rate,  even  though  the  predicted 
temperatures  also  averaged  1.3°C  colder,  which  would  tend  to  slow  the  melt  rate.  The 
actual  observed  iceberg  length  changes  as  compared  with  the  model  predictions  made  as 
part  of  this  test  were  inconclusive  due  to  the  time  constraints  (i.e.,  2  to  6  days.) 

The  various  studies  suggest  that  the  iceberg  deterioration  model  is  a  reasonable 
approximation  of  the  deterioration  process  when  observed  environmental  data  is  used, 
although  no  iceberg  has  been  observed  to  the  point  of  100%  melt.  The  1987  IIP  study 
(Hanson,  1987)  identified  significant  differences  between  the  FNMOC  data  and  observed 
data.  In  1988,  all  FNMOC  environmental  products  were  improved  and  the  new  data  was 
used  by  IIP  (Hanson,  1988).  The  new  values  for  SST,  sea  height,  and  sea  period  are 
reportedly  in  agreement  with  observed  values,  although  no  validation  experiments  have 
been  reported.  None  of  the  IIP  analyses  or  reports  indicates  that  a  complete  sensitivity 
analysis  of  the  deterioration  model  has  been  conducted. 


SENSITIVITY  ANALYSIS  OF  THE  ICEBERG  DETERIORATION  MODEL 
Analytic  Sensitivity  Analysis. 

Equations  (1)  -  (5)  (with  =  2)  describe  the  daily  melt  ra/e  of  icebergs. 
Combining  equations  (1)  -  (5),  the  daily  melt  rate  is 

V=  0.0012878xi^  +0.0301928X,  +0.008905  +  +1)/100  + 

0. 14490 Ix2\x,+1)/X5  (6) 

and  FC  is  defined  following  equation  (4).  The  first  order  sensitivity  of  the  daily  melt  rate 
to  changes  in  parameter  values  can  be  obtained  directly  from  equations  (1)  -  (5)  after 
differentiation  with  respect  to  the  parameters  of  interest. 

With  respect  to  sea  surface  temperature  (x,),  the  change  in  melt  rate  is  given  by 

—  =  0025756x,+.0101928+— +  0.144901-^  (7) 

a-,  100  Xj 
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Equation  (6)  indicates  that  the  melt  rate  is  quadratic  in  the  value  of  the  sea  surface 
temperature  and  that  the  rate  increases  with  increasing  temperature  (equation  (7)).  This 
indicates  that  errors  overestimating  the  SST  have  a  more  rapid  melting  effect  than  do 
corresponding  underestimating  errors. 

With  respect  to  wave  height  (x^),  the  change  in  melt  rate  is  given  by 
—  =  0.1159208^^^1^  (8) 


From  equation  (6),  the  melt  rate  varies  as  the  0.8  root  of  wave  height.  Equation  (8) 
indicates  that  the  change  in  the  melt  rate  decreases  exponentially  with  the  wave  height. 
Tliis  means  that  if  the  wave  height  is  overestimated,  the  melt  rate  is  also  overestimated 
with  the  percent  overestimation  increasing  at  a  decreasing  rate. 

With  respect  to  relative  speed  (Xj),  assuming  Xj  <  25,  the  change  in  melt  rate  is 
given  by 

^  =  %^(0«4  -  0.2021og,„l4)  (9) 


The  melt  rate  varies  linearly  with  respect  to  changes  in  relative  speed.  Equation  (9) 
indicates  that  this  rate  of  change  is  constant. 


With  respect  to  waterline  length  (X4),  the  change  in  melt  rate  is  given  by 


K. 

<^4 


(>Vl+l) 

100 


(-0.202^) 

X4 


(10) 


The  melt  rate  varies  negatively  with  respect  to  the  logjo  of  the  waterline  length.  The  rate 
of  change  of  the  melt  rate  varies  negatively  in  inverse  proportion  to  the  waterline  length. 
Aji  iceberg  with  a  longer  waterline  length  melts  slower  than  does  one  with  a  shorter 
length.  Overestimating  the  length  of  an  iceberg  keeps  it  in  the  system  longer  because  there 
is  more  to  melt  and  it  melts  slower  because  it  is  longer. 

Finally,  with  respect  to  the  wave  period  (Xj),  the  change  in  melt  rate  is  given  by 
dY  -0.144901x2'^(xi+1) 


Equation  (6)  indicates  that  melt  rate  varies  inversely  with  the  wave  period  and  equation 
(1)  indicates  that  the  change  in  the  rate  varies  negatively  in  inverse  proportion  to  the 
square  of  the  wave  period.  A  shorter  wave  period  results  in  a  faster  melt  (equation  (6)). 
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The  above  analytical  analysis  of  the  model  equations  provides  insights  into  the 
direction  of  the  changes.  In  order  to  assess  the  magnitude  of  changes  in  the  output  (daily 
melt  rate)  with  respect  to  changes  in  the  input  parameters,  it  is  necessary  to  conduct  an 
empirical  evaluation. 

Empirical  Sensitivity  Analysis. 

Methodology. 

The  model  was  initially  inspected  to  determine  the  factor  or  factors  believed  to 
constitute  the  major  contribution  to  deterioration.  An  initial  noininctl  range  sensitivity 
analysis  was  performed  to  determine  the  factor  that  was  dominant  in  the  majority  of 
different  scenarios.  Once  this  dominant  factor  was  discovered  its  value  alone  was  varied 
to  present  three  different  "nominal"  scenarios.  Witliin  these  scenarios,  further  univariate 
parametric  sensitivity  analysis  was  performed  to  determine  how  uncertainties  in  the  other 
factors  propagate  through  the  model.  This  information  was  captured  by  perturbing  their 
values  from  the  "nominal"  within  a  plausible  range.  The  relative  contributions  to  the 
uncertainty  in  the  output  are  then  determined. 

Results. 

The  factors  wliich  constitute  the  major  contribution  to  uncertainty  are  sea  surface 
temperature  x,(SST),  and  wave  height  (XAMP).  A  nominal  range  sensitivity  analysis 
was  performed  by  changing  their  values  by  ten  percent.  It  was  found  that  x,  was  the 
dominant  factor  in  the  majority  of  scenarios.  However,  it  should  be  noted  that  at  low 
values  of  x„  Xj  was  dominant.  Therefore,  three  different  SSTs  were  used  to  dictate  the 
three  different  "nominal"  scenarios:  1,  6,  and  15°C.  The  nominal  values  of  the  other 
factors  remained  constant  from  one  scenario  to  another: 

X2  =  XAMP  =  6  ft  (182.9  cm) 

X3  =  RELSPD  =  25  cnVs 
X,  =  RLEN  =  100  m 
x,  =IPER=  10  s 
Xg  =  ZTIME  =  2  half-days 
X,  =  WEATHER  =  1 

Only  X2,  X3,  and  Xj  were  perturbed  uniformly  by  fifty  percent,  forty  percent,  and 
fifty  percent  respectively,  while  x^  (waterline  length)  was  varied  by  classification  and  with 
much  more  disparity.  In  all  cases,  misclassification  obviously  had  the  greatest  effect  on 
uncertainty  in  melt  rate.  Therefore  a*2,  X3,  and  X5  were  considered  together  in  teims  of 
their  relative  uncertainties. 
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Results  show  that  at  a  sea  surface  temperature  of  1°C,  changes  in  wave  height  had 
the  greatest  effect  on  the  deterioration,  followed  by  wave  period.  In  fact,  the  uncertainty 
in  wave  height  propagated  rouglily  tliirty  percent  more  uncertainty  as  wave  period. 
Relative  speed  had  an  almost  negligible  eflfect.  Results  at  a  sea  surface  temperature  of  6°C 
were  similar.  Uncertainty  in  wave  height  propagated  roughly  twenty  percent  more 
uncertainty.  Again,  changes  in  relative  speed  had  an  almost  negligible  effect.  At  a  sea 
surface  temperature  of  15°C,  results  were  almost  identical  to  those  at  6°C.  Detailed 
results  are  included  in  Appendix  C  in  Armacost  (1994). 

Table  1  illustrates  the  impact  of  the  joint  variation  in  sea  surface  temperature  and 
wave  height  on  melt  rate.  The  nominal  values  are  indicated  in  boldface.  Alone,  a  10% 
variation  in  sea  surface  temperature  (at  1°C)  results  in  about  a  4%  variation  in  melt;  a 
10%  variation  in  wave  height  alone  (at  6ft)  results  in  about  6.5%  variation  in  melt. 
However,  their  joint  10%  variation  results  in  about  12%  variation  in  melt. 

Table  1.  Ten  percent  parametric  variation  in  SST  (1°C)  and  Wave  Height  (6ft). 


SST  (“C) 

XAMP(cm) 

MELT(cm/day) 

0.9 

164.6 

1.924 

1.0 

164.6 

2.024 

1.1 

164.6 

2.125 

0.9 

182.9 

2.068 

1.0 

182.9 

2.176 

1.1 

182.9 

2.284 

0.9 

201.2 

2.208 

1.0 

201.2 

2.324 

1.1 

201.2 

2.439 

Table  2  includes  the  results  when  the  nominal  value  of  the  sea  surface  temperature 
changes  to  6°C.  Now  the  single  variation  in  sea  surface  temperature  results  in  8% 
variation  in  melt  rate  (due  to  the  larger  absolute  variation  in  SST)  wliile  the  wave  height 
variation  remains  around  6.5%.  The  joint  variation  now  averages  16%. 

Finally,  sensitivity  with  respect  to  size  classification  reveals  the  greatest 
opportunity  for  propagating  uncertainty.  Table  3  illustrates  the  time  to  100%  melt  for 
three  levels  of  sea  surface  temperature  with  other  parameters  held  constant  at  their 
nominal  values. 

Table  3  makes  it  very  clear,  particularly  at  cold  temperatures,  that  misclassifying  an 
iceberg  as  smaller  than  it  actually  is  will  result  in  the  iceberg  existing  long  after  it  has  been 
removed  from  the  plot,  even  if  waiting  until  125%  or  even  150%  of  melt  before  removing 
the  iceberg  from  the  system. 
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Table  2.  Ten  percent  parametric  variation  in  SST  (6®C)  and  Wave  Height  (6ft). 


SST  (°C) 

XAMP(cm) 

MELT(cm/day) 

5.4 

164.6 

6.470 

6 

164.6 

7.080 

6.6 

164.6 

7.691 

5.4 

182.9 

6.954 

6 

182.9 

7.609 

6.6 

182.9 

8.266 

5.4 

201.2  ^ 

7.428 

6 

201.2 

8.128 

6.6 

201.2 

8.829 

Table  3.  Classification  variation  and  time  to  100%  melt. 


Size 

SST  =  1“C 

SST  =  6°C 

SST  =  16“C 

Small  (60  m) 

26-27  days 

7-8  days 

3-4  days 

Medium  (122  m) 

55-56  days 

15-16  days 

6-7  days 

Large  (225  m) 

103-104  days 

29-30  days 

12-13  days 

Summary. 

The  empirical  analysis  confirms  the  importance  of  having  good  estimates  of  sea 
surface  temperature  and  wave  height.  The  results  confirm  what  has  been  known  about  the 
variation  in  melt  with  respect  to  changes  in  a  single  parameter.  A  new  result  from  tliis 
analysis  is  the  overall  variation  in  melt  with  respect  to  joint  variation  in  the  parameters. 
The  results  suggest  that  tliis  overall  variation  is  superlinear  (12-16%  output  variation  for  a 
10%  input  variation).  The  most  significant  result  is  the  impact  of  misclassification 

Simulation  analysis. 

The  above  parametric  analysis  provides  the  opportunity  to  isolate  sensitivity  effects 
with  respect  to  particular  parameters.  It  does  not  however  provide  an  ability  to  examine 
the  joint  effects  of  multiple  parameters  unless  specific  combinations  of  changes  are 
examined.  Clearly,  tliis  becomes  computationally  prohibitive  and  there  is  no  effective 
means  of  evaluating  the  resulting  outcomes.  An  alternative  means  of  examining  these 
effects  is  to  use  a  simulation  model  that  considers  the  parameters  to  be  random  variables 
with  specified  probability  distributions.  A  Monte  Carlo  simulation  then  can  determine  the 
distribution  of  an  output  variable  of  interest.  Unfortunately,  such  a  simulation  is  only 
descriptive  and  simply  describes  the  system  output  for  a  given  set  of  inputs.  It  does 
provide  the  capability  to  examine  various  inputs  of  interest  and  determine  how  the  system 
outputs  will  change. 

To  illustrate  the  capability  and  limitations  of  simulation,  the  deterioration  model 
was  simulated  in  the  Excel  environment  to  illustrate  the  use  of  simulation  to  evaluate  the 


Analysis  of  the  IIP  Iceberg  Deterioration  Model 


Page  9 


sensitivity  of  policy  variables.  Details  of  the  simulation  are  included  in  Appendix  I.  It  was 
assumed  that  the  input  variables  were  independent  and  normally  distributed  with  the 
parameters  indicated  in  Table  4.  Each  simulation  run  involved  28  half-days.  The 
simulation  involved  100  runs  (total  of  2800  half  days.) 

Table  4.  Simulation  parameters. 


Parameter 

Mean 

Standard 

Deviation 

Wave  Height 

182.9  cm 

30  cm 

Sea  surface  temperature 

6°C 

2°C 

Wave  period 

10  sec 

1 .66  sec 

Relative  speed 

25  cm/sec 

3.33  cm/sec 

The  purpose  of  the  simulation  is  to  examine  the  distribution  of  waterline  lengths 
after  28  half  days  in  order  to  evaluate  berg  deletion  policies.  The  28  half  days  corresponds 
to  the  approximate  revisit  cycle  of  the  IIP.  Assume  that  the  iceberg  drift  is  such  that  it 
remains  in  the  vicinity  of  the  LAKI  (60  nm)  during  that  period.  After  100  runs  for  a  small 
iceberg  (60  m  initial  waterline  length),  the  average  waterline  length  is  -45.17  m  with  a 
standard  deviation  of  7.46  m.  The  distribution  of  waterline  length  after  14  days  is 
illustrated  in  Figure  2.  The  distribution  is  approximately  normal. 


]  Frequency 
-  Cumulative  % 


Waterline  length 


Figure  2.  Waterline  length  for  a  small  iceberg  after  14  days  melt. 

With  an  initial  waterline  length  of  60  m,  a  150%  melt  deletion  policy  would  set  a 
waterline  length  of  -30  m  as  the  deletion  tlireshold.  Under  this  criterion,  97.9%  of  the 
small  icebergs  would  be  deleted  by  the  model  between  ICERECDET  patrols. 


Similar  simulations  were  conducted  for  medium  and  large  bergs.  The  resulting 
distributions  of  waterline  length  were  also  approximately  normally  distributed.  For 
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medium  bergs  (initial  waterline  length  of  122  m),  the  average  waterline  length  after  28  h^f 
days  was  12.31  m  with  a  standard  deviation  of  7.07  m.  Tliis  means  that  none  of  these  will 
have  been  deleted  by  the  model  (waterline  length  of -61m)  by  the  time  of  the  next  patrol, 
but  that  most  of  them  will  be  very  difficult  to  detect.  If  resighted,  they  should  be  classified 
as  a  smaU  berg  and  would  be  deleted  in  the  next  14  day  period.  The  probabUity  that  a 
medium  iceberg  wiU  remain  a  medium  iceberg  at  the  next  14  day  sighting  is  less  than  10“-\ 

For  large  bergs  (initial  waterline  lengtli  of  225  m),  the  average  waterline  length 
after  28  half  days  was  118.1  m  with  a  standard  deviation  of  6.38  m.  As  with  medium 
icebergs,  none  of  these  wiU  have  been  deleted  by  the  model  (waterline  length  of -122m)  by 
the  time  of  the  next  patrol.  Approximately  23%  will  still  be  classified  as  large  icebergs  and 
the  remaining  77%  will  be  classified  as  medium  icebergs. 

The  above  results  obviously  depend  on  the  nominal  values  of  the  input  parameters. 
Clearly,  these  change  over  the  IIP  area  of  responsibility.  In  one  area  a  nominal  SST  of  6°C 
is  reasonable  (virtually  all  observations  in  the  0-12°C  range  and  two-thirds  in  the  4-8°C 
range).  In  other  areas,  different  values  should  be  used.  Nonetheless,  the  vai  lability 
represented  should  more  than  adequately  capture  measurement  uncertainty.  The  above 
analysis  suggests  that  the  150%  deletion  policy  provides  good  protection  against  deleting 
an  iceberg  prematurely  in  those  areas  where  the  environmental  parameters  hold  and  the 
iceberg  drift  is  such  that  it  remains  in  the  area. 


SUMMARY  AND  CONCLUSIONS 

This  evaluation  and  analysis  of  the  IIP  iceberg  deterioration  model  has  concluded 
that  the  model  appears  to  be  a  very  reasonable  representation  of  the  deterioration  process, 
with  the  exception  of  calving  for  which  meaningful  data  will  be  impossible  to  obtain.  This 
conclusion  is  supported  by  our  review  and  that  of  others  described  previously.  The 
analytical  and  empirical  sensitivity  analyses  indicated  that  sea  surface  temperature  and 
wave  height  are  very  important  parameters  in  the  model  with  respect  to  their  effect  on  the 
output  (melt  rate,  time  to  melt).  However,  any  adverse  impacts  that  errors  in  these 
parameters  may  have  are  completely  overshadowed  by  the  effects  of  misclassification  of 
the  iceberg  (wrong  specification  of  the  initial  waterline  length).  The  simulation  analysis 
indicated  that  uncertainty  in  parameter  values  (including  sea  surface  temperature  and  wave 
height)  are  absorbed  by  the  deletion  policy,  assuming  that  the  iceberg  is  correctly 
classified.  Therefore,  it  does  not  appear  that  further  refinement  of  the  input  environmental 
variables  is  required.  Any  additional  effort  should  be  directed  toward  ensuring  a  correct 
initial  classification  of  the  icebergs. 
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Appendix  1.  Simulation  Evaluation  of  Iceberg  Deletion  Policies. 


This  Appendix  includes  an  analysis  of  iceberg  deletion  policies  using  a  Monte  Carlo 
simulation  of  the  IIP  iceberg  deterioration  model.  The  analysis  was  conducted  by  Richard 
Ashley,  Grisselle  Centeno,  Stephen  Joseph,  and  Lou  Tozer  at  the  University  of  Central 
Florida,  December,  1994.  The  value  of  this  analysis  is  as  a  demonstration  of  the  use  of 
simulation  to  conduct  sensitivity  analyses  with  respect  to  policy  variables  of  interest.  The 
conclusions  drawn  in  the  report  are  limited  by  the  assumptions  stated  therein. 
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EXECUTIVE  SUMMARY 


This  report  gives  an  overview  of  the  International  Ice  Patrol  (IIP),  including  its  origin, 
mission,  responsibilities,  and  issuance  of  the  Ice  Bulletins  and  Ice  Charts.  Also,  a 
description  of  the  flow  and  determination  of  information  necessary  for  the  execution  of  the 
deterioration  and  drift  models  is  included.  An  analysis  is  done  on  the  uncertainty  created 
by  the  inputs  to  the  deterioration  model,  how  this  uncertainty  is  propagated  throughout  the 
model,  and  how  it  affects  the  rules  governing  the  deletion  of  icebergs  from  the  model. 
Finally,  conclusions  based  on  the  analysis  are  made. 


IV 


INTRODUCTION _ _ 

On  April  15,  1912,  the  R.M.S.  Titanic  liit  an  iceberg  and^ank  in  the  trans-Atlantic 
shipping  lanes  near  the  Grand  Banks  off  Newfoundland.  On  January  20, 1914  the 
International  Convention  for  the  Safety  of  Life  at  Sea  was  signed  to  provide  for  the 
establishment  of  an  International  Service  of  Ice  Observation  and  Ice  Patrol.  Tire  current 
ice  patrol  authority  is  the  International  Convention  for  the  Safety  of  Life  at  Sea  (SOLAS). 
In  1974,  SOLAS  formed  the  International  Ice  Patrol  (IIP)  to  track  icebergs  and  provide 
warnings  to  mariners  of  the  potential  dangers  caused  by  the  presence  of  these  icebergs. 

The  mission  of  the  IIP  is  to  "determine  the  Limits  of  All  Known  Ice  along  the  southeastern, 
southern,  and  southwestern  edge  of  the  ice  region  and  publish  that  information  to  mariners 
in  a  timely  fashion". [1]  The  United  States  Coast  Guard  was  tasked  by  SOLAS  with  the 
management  and  operation  of  the  IIP.  The  area  of  the  HP's  responsibility  is  from  40T4  to 
52"N  latitude  and  39^  to  5TW  longitude.  Figure  I  shows  the  area  bordered  to  the  west 
by  Newfoundland,  Canada  and  to  the  northeast  by  Greenland. 

Over  10,000  icebergs  are  calved  from  the  western  coast  of  Greenland  each  year.  They  are 
then  caught  in  the  Labrador  Current  and  carried  southward  through  the  Davis  Strait  and 
the  Labrador  Sea  to  the  Grand  Banks  off  the  southeast  comer  of  Newfoundland.  Flowing 
north  then  eastward  off  the  northeast  coast  of  the  U.S.  is  the  much  warmer  Gulf  Steam. 
When  these  two  water  currents  collide  near  the  Grand  Banks,  their  contrasting 
temperatures  cause  dense  fog  to  form  nearly  half  of  the  year.  The  IIP  aggressively  tracks 
the  icebergs  that  cross  south  of  the  48"N  and  into  the  heavily  traveled  shipping  lanes 
during  the  iceberg  season  from  March  to  August. 

One  of  the  primary  products  of  the  HP’s  efforts  is  the  issuance  of  two  daily  Ice  Bulletins. 
The  OOOOZ  and  1200Z  Ice  Bulletins  are  radio  transmissions  disseminated  at  midnight  and 
noon  Greenwich  Mean  Time,  respectively.  NMF/NIK  the  radio  station  of  the  U.S.  Coast 
Guard  Communications  Station  in  Boston,  MA.,  and  VON,  the  radio  station  of  the 
Canadian  Coast  Guard  Radio  Station  at  St.  John's,  Newfoundland,  are  the  primary  radio 
stations  responsible  for  these  transmissions.  The  Bulletins  contain  the  following 
information:  the  estimate  Limits  of  All  Known  Ice,  the  estimated  limit  of  sea  ice,  positions 
of  southern  and  eastern  most  icebergs,  position  of  growlers,  positions  oj  radar  targets,  and 
the  area  of  many  icebergs.  Figure  2  shows  an  example  of  a  1200Z  Ice  Bulletin  transmitted 
July  21,  1994.  [1] 

The  other  primary  product  is  the  HP's  issuance  of  the  1200Z  Facsimile  Chart  which 
depicts  the  Limits  of  All  Known  Ice.  The  U.S.  Coast  Guard  Communication  Station 
NMF/NIK  transmits  the  ice  limits  chart  at  1600Z  and  1810Z  (or  4P.M.  and  6P.M. 
Greenwich  Mean  Time)  daily  during  the  iceberg  season  along  with  safety  messages  to 
warn  mariners  of  icebergs  sighted  outside  of  the  published  ice  limits.  Figure  3  shows  a 
1200Z  Facsimile  Chart  transmitted  July  21,  1994.  [1] 
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NC  DE  U9 

O  202250Z  JUL  94 
FM  COMINTICEPAT  GROTON  CT 
TO  AIG  EIGHT  NINE  ONE  SIX 

CANADIAN  COASTAL  RADIO  STATION  ST  JOHNS  NFLD 
MV  PAUL  BUCK 
COMSAT  SAFETYNET 

950  L  ENFANT  PLAZA 
WASHINGTON  DC  20024 
XMT  HMCS  HALIFAX 
BT 

UNCLAS  //N16170// 

SUBJ:  INTERNATIONAL  ICE  PATROL  (IIP)  BULLETIN 
SECURITE 

CCODE/1 : 31 ; 04 : 01+1200 : OO/AOW/ IIP/CCODE 

1.  211200Z  JUL  94  INTERNATIONAL  ICE  PATROL  (IIP)  BULLETIN. 

REPORT  POSITION  AND  TIME  OF  ALL  ICE  SIGHTED  TO  COMINTICEPAT  VIA  CG 
COMMUNICATIONS  STATION  NMF,  NMN,  INMARSAT  CODE  42,  AND  ANY  CANADIAN 
COAST  GUARD  RADIO  STATION.  ALL  SHIPS  ARE  REQUESTED  TO  MAKE 
UNCLASSIFIED  SEA- SURFACE  TEMPERATURE  AND  WEATHER  REPORTS  TO 
COMINTICEPAT  EVERY  SIX  HOURS  WHEN  WITHIN  THE  LATITUDES  OF  40N  AND 
52N  AND  LONGITUDES  39W  AND  57W.  IT  IS  NOT  NECESSARY  TO  MAKE  THESE 
REPORTS  IF  A  ROUTINE  WEATHER  REPORT  IS  MADE  TO  METEO  WASHINGTON  DC. 
ALL  MARINERS  ARE  URGED  TO  USE  EXTREME  CAUTION  WHEN  TRANSITING 

NEAR  THE  GRAND  BANKS  SINCE  ICE  MAY  BE  IN  THE  AREA. 

2.  THE  ICEBERG,  GROWLER,  AND  RADAR  TARGET  POSITIONS  ARE  BASED  ON 
ESTIMATED  DRIFT.  DATE  OF  SIGHTING  IS  IN  PARENTHESIS  FOLLOWING  THE 
POSITION.  ALL  DATES  ARE  JULY  UNLESS  OTHERWISE  INDICATED. 

3.  ESTIMATED  LIMIT  OF  ALL  KNOWN  ICE:  FROM  THE  NEWFOUNDLAND  COAST 
NEAR  4635N  5310W  TO  4425N  5100W  TO  4415N  4700W  TO  4500N  4300W  TO 
5100N  4035W  TO  5200N  4100W  TO  5700N  5500W  THEN  EASTWARD. 

THE  ICEBERG  LIMIT  NORTH  OF  52N  IS  OBTAINED  FROM  ENVIRONMENT 
CANADA'S  ICE  CENTER  OTTAWA. 

4.  NO  SEA  ICE  SOUTH  OF  52N. 

5.  SOUTHERN  AND  EASTERN  MOST  BERGS  ESTIMATED  AT:  4550N  4338W(16), 
4459N  4738W(15),  4450N  4743W(15),  4450N  4744W(15),  4536N  4358W(12), 
5114N  4134W(12),  4633N  4313W(06),  4454N  4817W(17). 

6.  RADAR  TARGETS  ESTIMATED  AT:  4721N  4358W(18),  4718N  4356W(18). 

7.  THE  FOLLOWING  RADAR  TARGETS  ARE  OUTSIDE  THE  LIMITS  OF 

ALL  KNOWN  ICE:  4627N  5415W(19),  4609N  5330W(19),  4705N  3830W(14). 

8.  THERE  ARE  MANY  ICEBERGS  AND  GROWLERS  NORTH  OF  4500N  AND  WEST 
OF  4400W  WITHIN  THE  LIMITS  OF  ALL  KNOWN  ICE. 

BT 


Figure  2:  2I1200ZJUL  94  International  Ice  Patrol  (IIP)  Bulletin 
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Figure  3:  1200  GMT  21  JUL  94  International  Ice  Patrol  (UP)  Facsimile  Chart 


LAKI  PROCESS  DESCRIPTION 


The  IIP  determines  the  LAKI  using  many  different  sources  of  information  in  the  generation 
of  the  Ice  Bulletins  and  the  Ice  Chart.  See  Figure  4  for  the,  information  flow.  A 
description  of  tire  input  information  follows. 

Size  and  shape  classification 

m  Over  60%  of  the  data  is  radar  based. 

■  Less  than  40%  are  visual  sightings  received  from  ships  operating  in  the  area. 

■  Four  categories  for  size:  growler,  small,  medium  and  large  (die  iceberg  is  assumed 
to  be  medium  when  it  is  considered  a  target  iceberg). 

■  Two  categories  for  shape:  tabular,  non  tabular  (the  iceberg  is  assumed  to  be  non- 
tabular  when  it  is  considered  a  target  iceberg). 

Environmental  information 

■  This  information  is  received  primarily  from  U.S.  Navy  Fleet  Numerical 
Meteorology  and  Oceanography  Center  (FNMOC  ). 

■  The  information  includes:  surface  wind,  wave  height,  wave  period,  sea  surface 
temperature. 

Water  current  information 

■  This  information  includes:  Geostrophic  mean  historical  current,  and  real-time 
(local)  received  from  drift  buoys  tracked  by  satellite. 

200%  Side  Looking  Airborne  Radar  (SLAR) 

■  The  Ice  Reconnaissance  Detachment  (ICERECDET)  conducts  a  200% 
surveillance  of  the  area  every  14  days.  This  information  is  used  to  update  the 
positions  of  the  known  icebergs,  to  add  any  icebergs  not  already  detected,  and  as  a 
criteria  for  the  deletion  of  an  iceberg. 

The  IIP  Information  Collection  Center  (ICC)  uses  these  as  inputs  to  the  Drift  and 
Deterioration  Models  described  below. 

Drift  model 

■  Inputs  used:  local  wind,  position,  size  and  shape,  geostrophic  current,  local 
current. 

■  Description:  Tlie  drift  model  updates  the  position  of  the  icebergs  as  they  move 
through  the  area  of  die  IIP  responsibility. 

Deterioration  model 

m  Inputs  used:  position,  size,  sea  surface  temperature,  wave  height,  wave  period. 
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■  Ways  of  deterioration:  insulation  (sun  heating),  buoyant  convection  (vertical 
circulation  of  tlie  water),  wind  force  convection  (drift  movement  tlirough  tlie 
water),  wave  induced  (wave  washing  of  the  subaerial  surface). 

■  Deletion  criteria:  The  iceberg  is  deleted  from  the  database  of  all  known  icebergs  if 
125%  of  the  original  waterline  length  has  melted.  If  tlie  iceberg  is  a  limit  setting 
iceberg  then  it  is  retained  in  the  model  until  150%  of  its  waterline  length  has 
melted. 

■  Description:  Predicts  the  degree  of  iceberg  linear  surfece  length  melt. 

IIP  produce  the  Ice  Bulletins  and  Ice  Chart  only  after  the  Deterioration  and  Drift  Models 
output  is  sent  back  to  the  ICC  and  combined  with  all  inputs. 

Final  Outputs 

m  Ice  Bulletins:  See  Figure  2  for  an  example  of  an  Ice  Bulletin. 

■  Ice  Chart:  See  Figure  3  for  an  example  of  an  Ice  Chart. 


Figure  5:  Project  Information  Flow 
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PROJECT  OBJECTIVE 


The  main  objective  of  the  project  is  to  verify,  within  99%  confidence,  that  150%  melt 
criteria  for  the  limit  setting  icebergs  is  a  reliable  policy  decision  for  deleting  an  iceberg 
from  the  overall  iceberg  database.  Figure  4  shows  the  overall  information  flow  for  the 
generation  of  the  LAKI  Ice  Bulletins  and  Ice  Charts.  The  dashed  portion,  detailed  in 
Figure  5,  represents  the  information  flow  of  the  area  where  this  objective  is  applied. 


PROJECT  ASSUMPTIONS _ 

The  following  assumptions  were  drawn  in  the  analysis  of  the  deletion  criteria. 

■  The  icebergs  are  always  melting. 

■  A  limit  setting  iceberg  is  defined  as  an  iceberg  within  60nm  north  and  east  of  the 
LAKI,  or  any  location  south  and  west  of  the  LAKI.  The  project  only  considers 
these  icebergs  which  are  consider  to  stay  within  this  limit  setting  area  throughout 
the  analysis. 

■  The  water  length  used  in  tlie  Monte  Carlo  simulation  is  the  maximum  waterline 
length  for  the  small,  medium,  and  large  size  icebergs. 

■  It  is  assumed  that  the  200%  SLAR  coverage  at  the  end  of  the  28th  half  day  did  not 
see  the  simulated  icebergs.  This  assumption  along  with  the  use  of  only  limit 
setting  icebergs  is  important  because  if  an  existing  limit  setting  iceberg  is  missed 
by  the  reconnaissance  and  deleted,  the  result  could  be  catastrophic. 

■  A  Monte  Carlo  simulation  is  run  for  28  half  days  (14  days),  this  time  is  considered 
discrete  time 

■  All  input  variables  for  the  Monte  Carlo  simulation  model  are  independent  and 
normally  distributed,  except  for  the  size  of  the  iceberg. 


UNCERTAINTY _ 

There  is  uncertainty  involved  at  all  levels  of  the  LAKI  process.  A  considerable  amount  of 
information  that  serves  as  inputs  for  the  drift  and  deterioration  models  are  derived  from  a 
historical  database.  This  along  with  measuremental  errors  contributes  to  the  uncertainty 
associated  with  the  models'  input  parameters. 

Otlier  sources  of  uncertainty  are  associated  with  the  initial  classifications  of  the  icebergs. 
For  example,  if  the  size  is  not  known,  then  it  is  assumed  to  be  medium.  As  a  result  of 
uncertainties  in  the  inputs,  the  model  outputs  are  also  uncertain.  The  propagation  of 
uncertainty  per  halfday  associated  with  the  deterioration  model  is  shown  in  Figure  6. 

There  is  also  uncertainty  associated  witli  position,  which  is  both  an  input,  and  an  output  of 
tlie  drift  model. 
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In  order  to  compensate  for  positional  uncertainty  the  IIP  uses  error  circles  which  increase 
in  diameter  over  time.  This  rule  also  takes  into  account  the  iceberg  splitting.  The  IIP  also 
uses  tlie  125%  and  150%  melt  criteria  for  deletion  which  means  that  an  iceberg  must  melt 
to  a  negative  waterline  length  before  it  is  deleted  from  the  all  known  iceberg  database. 

These  uncertainties  contribute  to  the  risk  of  collisions  between  ships  and  icebergs 
potentially  causing  human  and  monetary  losses.  Therefore  they  must  be  taken  under 
serious  consideration. 


METHODOLOGY _ _ _ _ _ 

Uncertainties  in  input  parameters  were  propagated  through  the  deterioration  model  using  a 
Monte  Carlo  Simulation.  This  method  was  chosen  because  imprecisions  in  the  propagated 
output  distribution  and  measures  of  uncertainty  are  easily  estunated.  There  is  also  no  need 
to  discretize  the  distributions  as  required  by  other  methods.  Improvement  can  be  achieved 
by  taking  more  samples  [5]. 

Once  a  preliminary  univariate  parametric  sensitivity  and  uncertainty  analysis  was 
performed  [2],  the  fectors  contributing  most  to  output  uncertainty  were  identified.  All 
input  distributions  are  considered  normal.  All  input  variables  are  randomly  generated  and 
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arc  assumed  independent.  The  base  case  values  were  used  as  expected  values  for  the  input 
variables.  The  plausible  uncertainties  in  the  input  variables  were  used  to  determine  die 
standard  deviations. 

The  predicted  melt  for  waterline  lengths  of  60m,  122m,  and  225m,  was  computed  using 
normally  distributed  random  variables.  This  melt  was  subtracted  from  the  previous  length 
and  used  as  the  new  length.  Then  another  set  of  random  variables  were  generated  to 
compute  die  next  half  day's  melt.  This  allows  for  the  combination  of  errors  over  all 
iterations.  Tliis  was  repeated  for  a  total  of  28  intervals  corresponding  to  14  days.  These 
iterations  constituted  one  simulation  run  (see  Appendix  II  for  spreadsheet,  and  Appendbc 
III  for  simulation  code). 

A  total  of  100  simulation  runs  were  performed  for  each  size.  Tlie  melt  was  combined  into 
intervals  for  a  histogram  to  determine  the  distribution.  This  was  then  verified  widi  the 
goodness  of  fit  using  the  chi-square.  The  uncertainty  in  the  output  distribution  represents 
the  uncertainty  associated  with  iceberg  melt  for  specific  inidal  lengths.  The  quantification 
of  these  errors  allows  for  the  determination  of  a  value  for  the  percent  melt  that  should  be 
used  as  a  criterion  for  iceberg  deletion. 


RESULTS 


The  results  for  the  100  runs  for  die  small,  medium  and  large  icebergs  can  be  found  in 
Appendix  I.  For  each  size  iceberg  descriptive  statistics  were  generated,  they  are  shown  in 
Tables  1,  5,  and  9  .  A  histogram  was  plotted  to  help  identify  the  proper  modeling 
distribution  (See  Figures  7,  8,  and  9).  For  all  cases  a  normal  distribution  appeared  to  be 
appropriate.  Based  on  the  normal  distribution,  the  distribution  parameters  from  the 
descriptive  statistics  were  used  to  develop  dieoretical  frequencies.  This  was  compared 
with  the  actual  frequencies  to  determine  Chi-Squared  values  for  a  Goodness-of-Fit  test. 
These  results  can  be  found  in  Tables  2,  6,  and  10. 

A  99%  upper  confidence  interval  on  tlie  standard  deviation  was  generated  for  each  size 
iceberg.  The  results  can  be  found  in  Tables  3,  7,  and  11.  These  upper  confidence  limits 
were  used  to  calculate  the  uncertainty  in  melt  for  each  given  size.  To  allow  for  an  iceberg 
to  be  deleted,  the  iceberg  must  be  melted  to  a  size  of  zero  plus  the  uncertainty.  Finally, 
tills  is  shown  as  %  melt  in  Tables  4,  8,  and  12  for  several  probabilities  of  existence. 
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Frequency 


For  Size  60  Berg 


Descriptive  Statistics 

Mean 

-45.17 

Standard  Error 

0.746218 

Median 

-43.8747 

Standard  Deviation 

7.46218 

Sample  Variance 

55.68412 

Kurtosis 

-0.42814 

Skewness 

-0.42338 

Range 

33.59649 

Minimum 

-64.242 

Maximum 

-30.6455 

Sum 

-4517 

Count 

100 

Larpest(l) 

-30.6455 

Smallest(l) 

-64.242 

Table  1:  Small  Iceberg  Melted  Size  Statistics 


Histogram 


Figure  7:  Small  Iceberg  Melted  Size  Histogram 
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For  Size  60  Berg 


Bin 

Frequency 

Cumulative  % 

Theor. 

Cumulative  % 

Theor. 

Frequency 

Chi- 

Square 

-64 

1 

1.00% 

.58% 

0.58 

0.30 

-60 

3 

4.00% 

2.34% 

1.76 

0.87 

-56 

6 

10.00% 

7.33% 

4.99 

0.20 

-52 

10 

20.00% 

18.00% 

10.67 

0.04 

-48 

12 

32.00% 

35.23% 

17.22 

1.58 

-44 

18 

50.00% 

56.23% 

21.00 

0.43 

-40 

22 

72.00% 

75.58% 

19.35 

0.36 

-36 

19 

91.00% 

89.04% 

13.46 

2.28 

-32 

7 

98.00% 

96.12% 

7.08 

0.00 

-28 

2 

100.00% 

98.93% 

2.81 

0.23 

-24 

0 

100.00% 

100.00% 

1.07 

1.07 

Total  -> 

7.37 

Chi-Square(.05,8)  ->  15.51 


Table  2:  Small  Iceberg  Goodness  of  Fit  with  Normal 


For  Size  60  Berg 


S* 

a 

Upper  C.i. 

S* 

Upper  C.L 
S 

55.68 

1.00% 

79.63 

8.92 

Table  3:  Small  Iceberg  99%  Confidence  Interval  on  Standard  Deviation 


For  Size  60  Berg 


P 

Existing 

Average 

Upper 

Limit 

Delta  Melt 
(error) 

%  melt 

.50% 

-45.17 

-22.18 

22.99 

138.31% 

1.00% 

-45.17 

-24.41 

20.76 

134.60% 

1.50% 

-45.17 

-25.81 

19.36 

132.27% 

2.00% 

-45.17 

-26.84 

18.33 

130.54% 

2.50% 

-45.17 

-27.68 

17.49 

129.15% 

3.00% 

-45.17 

-28.39 

16.78 

127.97% 

Table  4:  Probability  of  the  Existence  of  a  Small  Iceberg  v.s.  %  Melt 


Frequency 


For  Size  122  Berg 


Descriptive  Statistics 


Mean 

12.31291 

Standard  Error 

0.706614 

Median 

13.12535 

Standard  Deviation 

Sample  Variance 

7.066138 

49.93031 

Kurtosis 

-0.53565 

Skewness 

-0.13874 

Range 

31.6613 

Minimum 

-2.89096 

Maximum 

28.77034 

Sum 

1231.291 

Count 

100 

Largestd) 

28.77034 

Smallestd) 

-2.89096 

Table  5:  Medium  Iceberg  Melted  Size  Statistics 


Histogram 
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Figure  8:  Medium  Iceberg  Melted  Size  Histogram 


For  Size  122  Berg 


Bin 

Frequency 

Cumulative  % 

Theor. 

Cumuiative  % 

Theor. 

Frequency 

Chi- 

Square 

-2 

3 

3.00% 

2.14% 

2.14 

0.35 

2 

4 

7.00% 

7.22% 

5.08 

0.23 

6 

16 

23.00% 

18.58% 

11.36 

1.89 

10 

12 

35.00% 

37.17% 

18.59 

2.34 

14 

23 

58.00% 

59.44% 

22.26 

0.02 

18 

18 

76.00% 

78.95% 

19.52 

0.12 

22 

15 

91.00% 

91.48% 

12.53 

0.49 

26 

8 

99.00% 

97.36% 

5.88 

0.76 

30 

1 

100.00% 

99.38% 

2.02 

0.52 

34 

0 

100.00% 

99.89% 

0.51 

0.51 

Total  -> 

7.22 

Chi-Square(.05,8)  ->  14.07 


Table  6:  Medium  Iceberg  Goodness  of  Fit  with  Normal 


For  Size  122  Berg 


s* 

a 

Upper  C.l. 
S» 

Upper  C.l. 
S 

49.93 

1 .00% 

71.40 

8.45 

Table  7:  Medium  Iceberg  99%  Confidence  Interval  on  Standard  Deviation 


For  Size  122  Berg 


P 

Existing 

Average 

Upper 

Limit 

Delta  Melt 
(error) 

%  melt 

.50% 

12.31 

34.08 

21.77 

117.84% 

1.00% 

12.31 

31.97 

19.66 

116.11% 

1.50% 

12.31 

30.65 

18.34 

115.03% 

2.00% 

12.31 

29.67 

17.35 

114.22% 

2.50% 

12.31 

28.87 

16.56 

113.58% 

3.00% 

12.31 

28.21 

15.89 

113.03% 

Table  8:  Probability  of  the  Existence  of  a  Medium  Iceberg  v.x  %  Melt 
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For  Size  225  Berg 


Descriptive  Statistics 

Mean 

118.1247 

Standard  Error 

0.638016 

Median 

118.1752 

Standard  Deviation 

6.380162 

Sample  Variance 

40.70647 

Kurtosis 

-0.34303 

Skewness 

-0.15544 

Range 

28.07099 

Minimum 

102.0398 

Maximum 

130.1108 

Sum 

11812JZ 

Count 

100 

Largest(l) 

130.110^ 

Smallest(l) 

102.0398 

Table  9:  Large  Iceberg  Melted  Size  Statistics 


Histogram 


Figure  9:  Large  Iceberg  Melted  Size  Histogram 
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For  Size  225  Berg 


Bin 

Frequency 

Cumulative  % 

Theor. 

Cumulative  % 

Theor. 

Frequency 

Chi- 

Square 

102 

0 

.00% 

.57% 

0.57 

0.57 

106 

4 

4.00% 

2.87% 

2.29 

1.27 

110 

4 

8.00% 

10.14% 

7.27 

1.47 

114 

21 

29.00% 

25.90% 

15.75 

1.75 

118 

19 

48.00% 

49.22% 

23.32 

0.80 

122 

24 

72.00% 

72.82% 

23.60 

0.01 

126 

17 

89.00% 

89.15% 

16.33 

0.03 

130 

9 

98.00% 

96.86% 

7.72 

0.21 

134 

2 

100.00% 

99.36% 

2.49 

0.10 

Total  -> 

6.21 

Chi-Square(.05,8)  ->  14.07 


Table  10:  Large  Iceberg  Goodness  of  Fit  with  Normal 


For  Size  225  Berg 


S* 

a 

Upper  C.L 
S* 

Upper  C.l. 
S 

40.71 

1 .00% 

58.21 

7.63 

Table  II:  Large  Iceberg  99%  Confidence  Interval  on  Standard  Deviation 


For  Size  225  Berg 


P 

Existing 

Average 

Upper 

Limit 

Delta  Melt 
(error) 

%  melt 

.50% 

118.12 

137.78 

19.65 

108.73% 

1 .00% 

118.12 

135.87 

17.75 

107.89% 

1.50% 

118.12 

134.68 

16.56 

107.36% 

2.00% 

118.12 

133.79 

15.67 

106.96% 

2.50% 

118.12 

133.08 

14.95 

106.65% 

3.00% 

118.12 

132.47 

14.35 

106.38% 

Table  12:  Probability  of  the  Existence  of  a  Large  Iceberg  v.s.  %>  Melt 


CONCLUSION 


It  has  been  assumed  tliat  the  icebergs  in  the  model  have  not  been  located  by  the  SLAR. 
Based  on  this  assumption  and  the  uncertainty  determined  by  the  simulation  model,  IIP  s 
150%  melt  policy  is  a  conservative  policy  when  run  for  28  half  days  for  any  size  iceberg. 
For  example,  to  insure  an  iceberg  is  melted,  it  must  melt  100%  plus  the  uncertainty.  The 
simulation  model  shows  that  this  would  result  in  a  138%  melt  for  small  icebergs  instead  of 
the  existing  policy  of  1 50%  melt. 

Current  policy  states  that  if  the  SLAR  resighting  is  a  target,  the  current  deterioration 
model  size  is  used.  This  results  in  the  further  propagation  of  tlie  uncertainty.  To  study 
this,  a  simulation  was  run  for  56  halfdays  on  medium  icebergs  with  the  results  indicating 
an  increase  of  45 .7%  on  the  standard  deviation.  If  there  are  a  large  number  of  medium 
icebergs  that  are  resighted  as  targets,  this  percent  change  could  significantly  influence  the 

policy. 

The  uncertainty  in  melt  for  each  size  is  approximately  the  same.  (At  a  probability  of 
0.5%,  the  small  iceberg  uncertainty  was  22.99  m,  21.77  m  for  medium,  and  19.65  m  for 
large  icebergs.)  Therefore,  it  is  believed  that  the  policy  of  using  150%  melt  for  each  size 
iceberg  is  extremely  conservative  especially  for  die  larger  icebergs. 

By  lowering  tlie  melt  percent,  the  mariners  cost  could  potentially  be  reduced  without 
increasing  the  potential  danger  of  iceberg  collision,  if  the  area  of  the  LAKI  could  be 
reduced.  This  may  be  a  good  area  for  further  consideration.  For  example,  a  policy  of 
138%  for  small,  117%  for  medium,  and  108%  for  large  icebergs. 
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APPENDIX  I 


Simulation  results  of  100  runs  for  small,  medium,  and  large  icebergs. 
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For  Size  60 


Start  Size 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


60 


New  Size 


-47.1 


-41.5 


-54.8 


-36.6 


-42.0 


-61.3 


-42.9 


-40.7 


-41.4 


-40.8 


-39.1 


-46.3 


-53.3 


-36.6 


-42.7 


-48.8 


-56.3 


-35.5 


-55.7 


-47.7 


-51.5 


-49.0 


-41.0 


-39.6 


-44,2 


-53.2 


-57.4 


-34.4 


-39.3 


-53.6 


-38.3 


-43.2 


-39.4 


-41.1 


-58.7 


-47.1 


-51.0 


-46.9 


-46.8 


-39,3 


-64.2 


-49.2 


-41.1 


-38.3 


-48.8 


-42.7 


-55.0 


-47.8 


-52.6 


-45.1 


%Melt 


178% 


169% 


191% 


161% 


170% 


202% 


171% 


168% 


169% 


168% 


165% 


177% 


189% 


161% 


171% 


181% 


194% 


159% 


193% 


180% 


186% 


182% 


168% 


166% 


174% 


189% 


196% 


157% 


166% 


189% 


164% 


172% 


166% 


168% 


198% 


179% 


185% 


178% 


178% 


166% 


207% 


182% 


168% 


164% 


181% 


171% 


192% 


180% 


188% 


175% 


60 

-46.0 

60 

-40.0 

60 

-51.1 

60 

-48.3 

60 

-40.6 

60 

-61.2 

60 

-46.4 

60 

-41.8 

60 

-56.6 

60 

-54.4 

60 

-38.9 

60 

-33.5 

60 

-43.7 

60 

-30.6 

60 

-33.7 

60 

-47.4 

60 

-37.8 

60 

-58.9 

60 

-48.8 

60 

-36.2 

60 

-34.4 

60 

-41.9 

60 

-42.8 

60 

-45.0 

60 

-42.5 

60 

-44.8 

60 

-34.8 

60 

-52.3 

60 

-39.7 

60 

-50.7 

60 

-44.0 

60 

-47.9 

60 

-30.9 

60 

-48.0 

60 

-49.2 

60 

-60.1 

60 

-40.2 

60 

-51.6 

60 

-58.3 

60 

-38.3 

60 

-38.4 

60 

-43.3 

60 

-39.2 

60 

-39.8 

60 

-45.1 

60 

-34,6 

60 

-53.5 

60 

-42.0 

60 

-38.1 

60 

-39.0 

For  Size  122  Berg 


For  Size  225  Berg 


APPENDIX  III 


Simulation  program. 


Simulate  Macro 

Macro  recorded  11/19/94  by  Richard  Ashley 


Sub  Simulate ( ) 

For  I  =  1  To  Cells (3,  ”0") 

Range  (”F6").  Value  =  Range  (Cells  ( I  +  4,  ”R"),  Cellsd  +  4,  "R")). Value  ‘Transfer  new  berg  size  to  start  of  day  1 
For  J  1  To  1  'Loop  for  the  one  14  day  period 

Worksheets (  MeltSim” ). Range ( “RandNums” )• Calculate  'Generate  new  random  deviates 

WorksheetsC'MeltSim") . Range ( "Temp" ) .Calculate  'Recalc  T 

For  Half Day  =  1  To  28  'half  days 

Worksheets ("MeltSim" ) .Cells (HalfDay  +  5,  " F" ). Calculate 

Worksheets{"MeltSim") .Range(Cells(HalfDay  +  5,  "J"),  Cells (Half Day  +  5,  "O" )). Calculate 
Next  HalfDay 

Worksheets{'‘MeltSim"}  .Cells(29  +  5,  "F"  )  .Calculate 
Range  { "F6"  )  .Value  *=  Range  ("F34") 

Next  J 

Range (Cel Is { I  +  4,  "S")^  Cells(I  +  4,  "S")). Value  =  Range ( "F34" ) .Value  'copy  new  size  back 
WorksheetsC'MeltSim")  .Range(Cells(HI  +  4,  "T"),  Cellsd  +  4,  "U" )).  Calculate 
Next  I 

End  Sub 
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